Transient thermoelectricity in a vibrating quantum dot in Kondo regime 
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We investigate the time evolution of the thermopower in a vibrating quantum dot suddenly 
shifted into the Kondo regime via a gate voltage by adopting the time-dependent non-crossing 
approximation and linear response Onsager relations. Behaviour of the instantaneous thermopower 
is studied for a range of temperatures both in zero and strong electron-phonon coupling. We argue 
that inverse of the saturation value of decay time of thermopower to its steady state value might be 
an alternative tool in determination of the Kondo temperature and the value of the electron-phonon 
coupling strength. 
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Investigation of time-dependent electron transport in 
single electron transistors has gained considerable trac- 
tion lately as a result of tremendous advances in the 
burgeoning field of nanotechnology and their perceived 
potential to replace MOSFET transistors^ one day. De- 
velopment of quantum computers* 2 - and single electron 
guns 3 - are expected to benefit from advances in detection 
of electrons in real time too^. 

Transient current ensuing after sudden shifting of the 
gate or bias voltage 5 '-'® displays different time scale o 9 ' 10 . 
For an asymmetrically coupled system, interference be- 
tween the Kondo resonance and the sharp features in 
the contacts' density of states gives rise to oscillations 
in the long time scaled. Modeling the contacts' density 
of states in a realistic fashion via ab initio calculations 
yielded accurate predictions about transient curren t 12 ' 13 . 

Measurement of thermopower (Seebeck coefficient) S 
can provide additional insight into transport experiments 
because its sign is a valuable tool to determine the align- 
ment of orbitals of the impurity with respect to the Fermi 
levels of the contacts. To this end, significant progress 
has been achieved in performing thermoelectric measure- 
ments in molecular junctional— . 

Theoretical studies focused on incorporating strong 
correlation effects into this prototype. Kondo effect, aris- 
ing from a hybridization between the net spin localized 
within the impurity and the electrons in the contacts, is a 
prime example of strong electronic correlations. Schemes 
employing Ng's ansataA^ and Wilson's numerical renor- 
malization groupA 9 - both found that the sign of the ther- 
mopower can be adjusted by tuning the energy level of 
the quantum dot in the presence of Kondo correlations. 
When the electrodes are ferromagnetic, it was found that 
the thermopower is suppressed at low temperatures in 
parallel configuration and asymmetrically coupled an- 
tiparallel configuration 2 ^. 

Nevertheless, all of the aforementioned studies only 
took into account the steady state behaviour of thermo- 
electric transport and there was very little understanding 



of the temporal evolution of thermopower until now. In- 
deed, a recent work made the first attempt to elucidate 
the time-dependent behaviour of the Seebeck coefficient 
for a noninteracting system 2 ^. In this paper, we will take 
a step forward and study the transient response of ther- 
mopower for an interacting quantum dot suddenly moved 
into the Kondo regime with a time dependent gate volt- 
age. 

We can describe this device with Holstein Hamiltonian. 
Its three pieces representing the contacts, the quantum 
dot and the tunneling between them can be written re- 
spectively as 



H(t) = H c + H D {t)+H T {t) 



(1) 



where 
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The operators dl(d a ) and cl (cfc ao -) with a=L,R cre- 
ate (annihilate) an electron of spin a within the dot and 
in the left(L) and right(R) contacts respectively. Vk a 
and [i a represent the hopping amplitudes and chemical 
potentials whereas a^(a) creates(annihilates) a phonon. 
A is the strength of the electron-phonon interaction and 
loq is the phonon frequency. Throughout this letter, we 
will be using atomic units where Ti = ks = e = 1. Un- 
der the assumption that the hopping matrix elements are 
equal and have no explicit time and energy dependence, 
coupling of the dot to the contacts can be parameter- 
ized as r(e) = Tp(e) where T is a constant defined by 
f = 27r| V(e/) | 2 and p{e) is the density of states of the 
contacts. In the following, we will take parabolic density 
of states and same bandwidth for both contacts. 
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We will be interested in strong coupling regime. This 
amounts to resonant tunneling which entails longer elec- 
tron lifetime and strong electron-phonon interaction. 
Tunneling electrons create and destroy a phonon cloud 
and this leads to polaron formation at the junction. 
Phonon mode is not coupled to the leads because it is 
unrealistic to allow polaron formation in the leads which 
are made of metal. This model has long been used to 
study inelastic electron transport in molecular junctions 
in steady stated. 

Upon applying the unitary Lang-Firsov canonical 
transformation in order to eliminate the electron-phonon 
coupling term when the electron-phonon coupling is suf- 
ficiently strong compared to the tunnel couplings—, the 
dot Hamiltonian becomes 

H D (t) = SH D (t)S^ 

a 

+Udt.dfd\di + ujoa^a. (3) 

In this Hamiltonian, the dot level and the Hubbard inter- 
action strengths are renormalized as Idot(t) = £dot(t) — 
(A 2 /wo) and U = U — (2A 2 /wo) respectively due to the 
canonical transformation. When we apply the slave bo- 
son transformation to this Hamiltonian in U — > oo limit, 
double occupancy of the dot level is prevented. How- 
ever, standard diagrammatic techniques are not applica- 
ble anymore. We can overcome this problem by rewriting 
the original electron operator on the dot in terms of slave 
boson and pseudofermion operators as 



d a {t) = b\t)U(t) 
4(*) = 

subject to the restriction 

Q = b^b+J2.flu = h 



(4) 



(5) 



which guarantees the single occupancy of the dot level. 
The slave boson Hamiltonian turns out to be 



H(t) 



■ ujQa'a 



kaa 



(0) 



where the tunnel coupling is renormalized and given by 



Vka(t) = Vka(t)exp 



r a 2 / 




2 N Ph- - 


^)] 



(7) 



Using the slave boson and pseudofermion decomposi- 
tion of the original fermion operators on the dot, the 
retarded Green function can be expressed asS 2 - 

G R (t,h) = -^(t-ti)[G« eWo (t,t 1 ) J B<(t 1 ,t)e^*- t ) 
+G< eudo (t,t 1 )B R (t 1 ,t)e^-^} (8) 



In this expression, the phase factor is defined by 

4>(h-t) = -0[JV(1 - e-*^* 1 -*)) + 

(N ph . + l)(l-e i "°^ "*))], (9) 



where g is defined as g = — 
Einstein distribution N p h. = 



and N p h., given by Bose- 
1 function, represents 



the average number of phonons for temperature T and 
phonon frequency ujq. 
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FIG. 1: This figure shows the density of states of both 
contacts and the quantum dot in the initial and final states 
schematically. The temperature gradient between the con- 
tacts is also depicted. 

Double time retarded and lesser Green's functions 
for pseudofermions and slave bosons are determined by 
solving coupled Dyson equations in a cartesian two- 
dimensional grid. We note that the phonon phase fac- 
tors remain attached to the pseudofermion Green's func- 
tions in Dyson equations so that the electron-phonon 
interaction is accounted for properly 2 ^. The only re- 
maining ingredient to obtain a closed set of equations 
is the self energies. We resort to non-crossing approx- 
imation(NCA) to express the pseudofermion and slave 
boson self-energie s 10 ' 24 . NCA is known to provide accu- 
rate results for dynamical quantities except temperatures 
below T/Tk ~0.1 or finite magnetic fields. We will avoid 
these regimes here. The values of the double time Green 
functions are kept in a square matrix which is propagated 
diagonally to describe the time evolution in an accurate 
way. 

In linear response, conductance of the device is given 
by 



G(t) = *f 
and the thermopower can be expressed as 

£ia(t) 



S(t) 



TL\i(t) ' 



(10) 



(11) 
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where the Onsager relations are 
L n {t) = Tx 

*"(/> I ' ^mo^-M)^ 

and 

L 12 (t)=T 2 x 

Eq. ([11]) alongside with Eq. ([H]) and Eq. (H3J) is the 

central result of this work and we will explore its con- 
sequences in the following. We previously studied the 
behaviour of instantaneous conductance G(t) for a vi- 
brating quantum dot suddenly shifted into the Kondo 
regime in detail^. It is our intention in this paper to 
extend this analysis and shed light on the instantaneous 
thermopower S(t) for the same system. However, we will 
need to stay in linear response since the Onsager relations 
are valid only in this regime. 

Kondo effect is a leitmotif of many body physics oc- 
curring as a result of hybridization between the net spin 
localized inside the quantum dot and the continuum elec- 
trons in the leads. When the dot level is situated below 
the Fermi level at low temperatures, Kondo resonance 
emerges as a very sharp resonance pinned to the Fermi 
level of the contacts. A low energy scale called Kondo 
temperature provides a good estimate for the linewidth 
of the Kondo resonance. Kondo temperature is denoted 
with Tk and given by 

— (f) 4 =p(-^). ,14) 

In Eq. (fT4"|) D is the half bandwidth of the conduction 
electrons while T — Tp(ef). 

We will consider the behaviour of the instantaneous 
value of the thermopower immediately after the dot level 
is switched from t\ = — 51" to e 2 = —2T at t=0 via a 
gate voltage. A transition is triggered from a non-Kondo 
state to a Kondo state as a result of this abrupt move- 
ment. The density of states of the dot both in initial and 
final levels as well as the parabolic structure of the den- 
sity of states of the contacts is shown schematically in 
FigfT] Since it takes considerable amount of time for the 
Kondo resonance to form in the final stated, dynamical 
quantities would require a similar amount of time to ad- 
just and reach their steady state values. Consequently, 
their evolution in the long timescale is non-trivial and 
requires a careful analysis. 

Instantaneous thermopower is displayed for various 
ambient temperatures in Fig. Rafter switching to the fi- 
nal dot level in infinitesimal bias for zero electron-phonon 
coupling. In this figure, we clearly observe that the ther- 
mopower starts decaying from zero to its steady state 
value for all temperatures studied. After a certain decay 
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FIG. 2: This figure shows the instantaneous thermopower 
S(t) immediately after the dot level has been moved to its final 
position for 7\=0.0035r (red solid), T 2 =0.0028r (blue short 
dashed), T 3 =0.002ir (orange dot dashed), T 4 =0.0014r (black 
dotted), T5=0.0007F (green long dashed) in linear response 
without any electron-phonon coupling. 



time, it reaches its steady state value. This steady state 
is negative for all temperature values. However, its abso- 
lute value increases until a certain temperature. Once we 
increase the temperature beyond that, the steady state 
value starts decreasing but the decay time stays constant. 
The temperature at which the steady state thermopower 
reaches its minimum value has been previously identi- 
fied as the Kondo temperature^ and our results are in 
agreement with this observation. The extra insight our 
results provide is the saturation of the decay time of ther- 
mopower and inverse of the saturated decay time. This 
non-trivial and counter-intuitive result is shown quite 
strikingly in Fig. 2) We will elaborate on its microscopic 
nature later on. Both the temperature below which satu- 
ration occurs and the inverse of the saturated decay time 
are equal to the Kondo temperature. These novel ob- 
servations can become a powerful tool to determine the 
Kondo temperature during experiments involving ther- 
moelectric switching of a single electron transistor. 

In Fig. [31 we investigate the instantaneous ther- 
mopower for the same parameters used in Fig. [21 but 
this time we use g=2.25 with wo=0.08r. In this case, 
the dot level is shifted slightly downwards due to renor- 
malization leading to a smaller Kondo temperature com- 
pared to g=0 case. As a result, steady state thermopower 
values are smaller than g=0 case if the ambient temper- 
ature is above both Kondo temperatures. However, they 
are larger than g—0 case when the ambient temperature 
is lowered below both Kondo temperatures. 

Kondo temperature can once again be identified as the 
temperature at which steady state thermopower reaches 
its minimum value. This is also the same temperature 
below which decay time saturates in analogy with g=0 
case. The major difference with g=0 case lies in the 
fact that the saturated decay time is now much longer 
as one can clearly see for the lowest curves in both 



4 



0.0 
-0.1 
-0.2 
m-0.3 
-0.4 
-0.5 



v, 
'■•\ 

%\ 

\\ ^- 

w 

100 200 300 400 500 600 700 

rt 



FIG. 3: This figure shows the instantaneous thermopower 
S(t) immediately after the dot level has been moved to its final 
position for Ti=0.0035r (red solid), T 2 =0.0028r (blue short 
dashed), T 3 =0.002ir (orange dot dashed), T 4 =0.0014r (black 
dotted), T 5 =0.0007r (green long dashed) in linear response 
for p=2.25. 



Fig.[3]and Fig.[3]because the Kondo temperature is much 
smaller as a result of renormalization of the dot level. In 
fact, our calculations show that the decay time of ther- 
mopower scales with 1/Tk- That is why curves with 
and without electron-phonon coupling perfectly overlap 
in Fig. [4] This scaling has been previously observed for 
time dependent conductance calculations in the absence 
of any electron-phonon coupling as well-. Same under- 
lying mechanism holds for thermopower as well because 
both dynamical quantities depend on the evolution of 
Kondo resonance whose development has been shown to 
be inversely proportional to the Kondo temperature^. 
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FIG. 4: This figure shows the decay time of thermopower r 
as a function of temperature for both zero and finite electron- 
phonon coupling strengths. Saturation of the decay time of 
thermopower around Tk is clearly visible here. 

Electron-phonon coupling strength A is generally un- 
known in an experiment. These type of measurements 
may pave the way for determination of A if the phonon 



frequency uiq is known. Once the Kondo temperature 
is identified via the saturation of decay time described 
above, Eq. (fT4"]) can be solved for renormalized dot level 
e. This would in turn enable to find A because the final 
bare dot level is adjusted via gate voltage thus known. 

It is possible to understand previous numerical results 
intuitively with the help of Sommerfeld expansion at low 
temperatures. It is given by 



S(T) 



ir 2 T dA 
~3A(0,T)~£k 



c=0 



(15) 



in atomic units where A(0, T) is the value of the spectral 
function at Fermi level and ^ is its derivative. Since 
the Kondo resonance lies slightly above the Fermi levels, 
the derivative of the dot density of states at Fermi level 
is positive. This results in negative thermopower values 
for all temperatures in Kondo regime due to the negative 
sign at the beginning of Sommerfeld expansion. 

When T< Tk, the Kondo resonance is fully formed in 
steady state and lowering the temperature any further 
doesn't alter its final shape. Similarly, evolution to that 
final shape takes same amount of time for any T below 
Tk- Hence, the decay time of thermopower to its steady 
state stays constant. However, the steady state value 
keeps going down in magnitude because of T prefactor 
in Sommerfeld expansion. This prefactor becomes the 
only variable below Tk in Sommerfeld expansion gov- 
erning the steady state value since the final shape of 
spectral function is the same. Subsequently, steady state 
thermopower curve is nonmonotonic and goes to zero at 
T=0 eventually. This issue has been confirmed by other 
sources too^. On the other hand, the Kondo resonance is 
simply nonexistent when the dot is in its initial level since 
T^> Tk- Consequently, dot density of states is essentially 
flat at Fermi level giving rise to zero slope which implies 
zero thermopower. This is why S(t)=0 at t=0 and this is 
in good agreement with previous steady state results^. 

In conclusion, we investigated the temporal evolution 
of the thermopower in response to an abrupt movement 
of the dot level to a position where the Kondo effect is 
present. We identified both the temperature below which 
decay time saturates and the inverse of the saturated de- 
cay time as new tools to identify the value of the Kondo 
temperature and the electron-phonon coupling strength. 
Transient dynamics of thermopower offers a complemen- 
tary picture to the previous steady state results because 
providing a detailed analysis of turning this device on 
and off is crucial for practical device characterization. 

In principle, it should be possible to reproduce short 
time behaviour of our results with time dependent den- 
sity matrix renormalization group method, however we 
do not believe it would be sufficient to capture the long 
timescale governing the development of Kondo resonance 
since the number of states required to keep the trunca- 
tion error fixed grows exponentially with time making 
calculations involving long times unfeasible^. 

We believe that it is possible to perform this exper- 
iment with present day technology because latest ad- 
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varices in ultrafast pump-probe techniques enable to mea- 
sure the transients down to femtosecond timescal o 27 ' 28 . 
The entire set-up shown in FigQ] is kept in a dilution 
refrigerator in order to access the ambient temperatures 
around Tr . This means that both electrodes are in prin- 
ciple at same ambient temperature. One needs to induce 
a temperature gradient between the contacts in this ex- 



periment by irradiating one of the contacts with a laser 
beam. Consequently, we aim to motivate further experi- 
ments in this field with this paper. 
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